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ABSTRACT 


: This report describes the development of two methods of deter- 
mining the on-axis microwave power density in both the Fresnel and 
Fraunhofer regions, for circular-aperture antennas with any arbi- 
trary, axially symmetric aperture illumination. Both methods have 
been programmed for digital computation. These methods are described 
and are applied to a number of different aperture distributions which 
could conceivably occur in practice. The resultant on-axis fields 
are compared and some conclusions relating the aperture distribution 
and the resulting on-axis fields -are made. 
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INTRODUCTION 


The microwave fields associated with the high power radar sys- 
tems of today,-and the even higher power radars expected in the 
future, present definite hazards to personnel, ordnance and fuels 
unless appropriate safety precautions are incorporated. It is, 
therefore, necessary to determine which regions of the various 
radar fields exceed established safe working levels of power den- 
sity. The hazardous areas should preferably be examined through 
actual measurements with the radar in operation. However, such 
measurements are time consuming, costly and susceptible to many 
variables. It thus becomes essential that an accurate and reliable 
method of computing the theoretical power density on the axis of the 
beam be available both to facilitate optimum test planning and to 
aid in interpretation of experimental results. 


In general, resort to a highly sophisticated theoretical 
approach is not necessary in order to estimate the maximum values - 
of power density which may be encountered along the beam axis, which 
is the usual requirement when safety measures are the primary con- 
sideration. For such purposes, it is permissible to base calcula- 
tions on the simplifying assumption that the antenna aperture is 
uniformly illuminated by the feed. This assumption yields calcu- 
lated values of power density in the Fraunhofer (far) field which 
are at least as high as produced by any other distribution of the 
aperture field, hence the error introduced, if any, is on the side 
of safety throughout the far field. Moreover, this assumption 
results in predicted values of the power density maxima in the 
Fresnel field which, except for a few unusual cases, are not signi- - 
ficantly lower than those based on other distributions. For these 
reasons, the calculation of "safe" distances for circular microwave 
antennas in connection with personnel hazards studies at the Naval 
Weapons Laboratory has been based on the assumption of a uniformly 
illuminated aperture. Under this assumption, far field computations 
are reduced to a straightforward application of the inverse square 
law, and the prediction of maximum hazards in the Fresnel region is 
made possible through formulas which approximate the upper envelope 
of the power density contour in this region. References (a) and (b) 
describe these methods in detail. 


While methods based on uniform illumination are useful for 
establishing safe distances and other cases when only the Fresnel 
field maxima are of interest, they are not sufficient for a compre~« 
hensive evaluation of measurements obtained in the Fresnel field of 
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non-uniformly illuminated antennas since instrumentation might be 
located at points of minimum or intermediate ‘intensity. While 
uniform illumination theoretically gives maximum antenna gain, 
practical antenna design usually necessitates reduction of the 
power at the aperture edge to decrease the loss from "spill-over" 
and decrease the power in the sidelobes. Thus the aperture ampli- 
tude distribution in practical antennas is usually designed so that 
the power is concentrated at the aperture center and is tapered to 
lower values at the edge. 


The effect of these tapered illuminations on the axial field 
has been calculated by other investigators for a few distributions 
expressible as analytic mathematical functions. . Not all distribu- 
tions of interest can be adequately approximated by simple func- 
tions, however. Even for those distributions which can be so 
expressed, the manual solution of the Fresnel equation for a suf- 
-ficient number of points is. time-consuming and difficult. It was 
thus desirable to develop a method by which the axial field could 
be conveniently calculated by computer for any arbitrary aperture 
distribution. Two somewhat similar methods, each having certain 
advantages, have been developed to fulfill this requirement. This 
report describes these methods and their application to the common 
distributions generally used to approximate the aperture illumina- 
tion of practical antennas. The methods have also been applied to 
a number of arbitrary distributions for which the on-axis field 
could not readily be calculated by the usual means. 


THEORY 


Consider the antenna aperture of Figure 1. The scalar dif- 
fraction field at any point P located at a distance R along the 
axis of a circular aperture of constant phase is given by the 
Fresnel approximation equation, (Reference (b), Page 198). 


ent a 
nee al ; F(p, 6) end" P/AR p dp ag! (1) 
fe) re) 
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where 


U, is the field at a point P on the axis of the aperture 


R is the axial distance from the aperture to the 
point P 


X is the wave length of the radiation 
@ is the radius of the aperture 
po is the radial coordinate in the plane of the aperture 


g’ is the angular coordinate of a point in the.plane of 
the aperture 


F(p, ¢’) is the field strength or amplitude taper in the plane 
of the aperture 


k .is the wave number on/2 


For the case of axially symmetric aperture illumination, 
equation (1) reduces to 


a : : 
-jkR 2 
i) F(p) eI" PAR p ap (2) 
Le} 


, ke 
Un = 9 


if we now define the dimensionless parameters 


we have 
: 2 
-4kR: . 
=j jee | F(r) e7d™ * [x r dr (3) 
x 
- . 


If the function F(r) is normalized with respect to its value at 
the center of the aperture, the magnitude of Up will be the ratio of 
the field at the point P to that at the centes of the aperture. 
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For the special case of uniform aperture illumination, 
F(r) = 1.0 and equation (3) reduces to 


1 
-jkR . 
U4 ee | endt 2°/x » ar (4) 
° 
= j e JER E cin] ~3n/ax (5) 
2x : 


Method of Finite Differences 

One approach to obtaining a digital computer solution for the 
case of axially symmetrical illumination of the aperture lies in 
converting equation (3) directly into a finite difference equation. 


Equation (3) may be rewritten to eliminate the exponential 
under the integral sign: : 


Up seegen[ f F(r) ny con (#2) 2 sf F(r) sin (2° =) x er]. (3a). 


By taking the absolute value of both sides of equation (3a) and squar-~ 
ing, we obtain an expression for the ratio of the power density Wp at 
the point on the aperture axis to the illumination power density Wo at 
the center of the aperture. Thus 


= lu, - (2 s) {Lf f cos (SE): x] 
-( ‘Y xz) ain(S2)z ae] } 


(6) 
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If we divide the aperture radius a into N equal increments 
Ae;, where Ary designates the ith increment from the periphery of. 
the aperture (see Figure 2), we have 


yr, <Nt12-4i ond at oe, -2 


N 


Further, letting e, represent the ordinate of the normalized 


aperture field intensity distribution at the center of the inter- 
val Mr,, we have 


to {LE yeas) ow seen] 


N 
(7) 
- ; / eae s 
+] 5 eo, (HtMe= 1) sing (He ue=2)| i 
a N a < N 


=1 


It is convenient to express the dimensionless power density 
ratio in terms of the average power density over the aperture W,, 
rather than in terms of Wo. 


Now 


(8) 


Substitution of this into equation (7) yields the following 
‘expression: 


2 
OO ath ni et On sn NOG (tue - 1) cos F(X U2 = 3) 
N a N 
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Method of Superposition 


An. alternative approach follows from the concept that the 
aperture amplitude distribution may be approximated by a number 
of uniform fields distributed over coaxial disks or "sub-apertures" 
of varying radius. At each point on the aperture axis, the solu- 
tions for each of the uniformly illuminated sub-apertures may then 
be combined, by the principle of superposition, to yield the net 
power density at the field point. This concept is illustrated in 
Figure 3. 


For each of the sub-apertures, equation (5) takes the form 


2 2 2 
mr wr . mr ; 
Ue). Se 2ein ss e6s a g sin S14) Le (10 
p’4 2x 2x 2 : 


where the subscript i pertains to the it® sub-aperture, and (U ot 
is the field at any point P on the aperture axis due to the 
ith sub-aperture, Ae; = (e1 - e4-1) is the amplitude of the field 


associated with the ith uniformly illuminated sub-aperture, and 
r,; is the radius of the ith sub-aperture. 


Note that the phase angle component j e~JKR nas been dropped. 
Since this angle is independent of r,, it will not alter the rela- 
tive phase of the contributions of the various sub-apertures and 
hence serves no useful purpose in this analysis. 


If the difference in radius between successive sub-apertures 


is chosen to be aconstant, = and if the i's progress from the 
N+i1-i 
N 


Substitution of the above quantity in equation (8) yields 


periphery toward the center, we have Yr, = 


2 2 
= _ ct f/N+1i-i x (N+1-i 
(U,) 4 =2 fe, sin = Raic2) | coe oe (43-4) 
2 
oe ioe fs cect gee 
vo sess)] 


(11) 
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METHOD OF FINITE DIFFERENCES 


NORMALIZED APERTURE RADIUS 


Figure 2 


Approximation of aperture distribution by finite increments with N = 10. 


NORMALIZED APERTURE AMPLITUDE A(r) 
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Summing the contributions of the N sub-apertures, 


2 2 


N N i . ; : ; 
zx (U,), =2| 2 te, sin (HA2= 4) dos (H22=4) 
i=l. P i=l ex N Ox N — 


. (12). 
N e- 2 
-j © fe, sin? £ Heit) | 
. T=]. -% ex N 
Taking the absolute value and squaring gives 
2 


W. N . i N 2 2 
B=! 5 (U_) =4 = de, sind Hei-4) con & (H#i= 3) 
i=1 2x N 2x N 


N 2 
+| 5 2 cin? & (Uti) . ie) 
i= x N 

In this method 
2 2 
N N 
We 2 5 ot g(Htiat) 2 F oe (Heid) (14) 
Wo ®an 7 N faq # N 
where 
Be oe ee : 
for Sy Spa 


Finally, upon substituting equation (14) into equation (13), 
we have . 


isl 


. 2 
ae i ru Ntles +l « eee 42-1) 
ze, Oy.) 8D oF _N 2x N ) 


. 2 : 
N : c 2 ; (15) 
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Equations (9) and (15) were programmed for the TEM. 7090 and 
STRETCH (IBM 7030) computers and represent the two methods of 
obtaining the on-axis fields of arbitrary aperture distributions. 


While the method of superposition was found to give more 
accurate results than the finite differences method for a given 
value of N, either method will give results as precise as desired 
by choosing a large enough value for N. It is desirable, however, 
to use as low a value of N as is consistent with the necessary 
accuracy, since each e, must be taken either from a graph of the 


aperture distribution (if it is irregular) or calculated from the 
mathematical formula for the distribution (if it can be represented 
by a simple mathematical function). In the following calculations 
a value of N = 50 was used. 


While the method of superposition is more accurate, it does not 
‘lend itself well (although it- could be modified to do so) to calcula- 
tions of the field when the aperture distribution is not single valued 
or monotonic, e.g. Taylor distributions, distributions with blocking, 
etc. For distributions of this type, the method of finite differences 
is especially useful. 


In this report the method of superposition was used to calculate 
the fields of all single valued distributions, with the method of 
finite differences being incorporated for the remaining aperture 
illuminations. In all the calculations it is assumed that the field 
is in phase across the entire aperture. 


DISCUSSION 


The aperture distributions for which the on-axis fields are 
calculated are uniform illumination; linear, square law, (1 - r=)? 
and (cosine)" tapers; Taylor distributions; some irregular arbitrary 
distributions and distributions with "blocking", where the field 
strength or power density is zero near the center of the aperture to 
simulate the effect of the antenna primary feed blocking the center 


“illumination. For the linear, square law ((1 = y?)" taper with n = 1) 
and (cosine)" tapers, four variations of each taper are used, i.e., 
tapers in which the edge illumination is 5, 10, 20 and » db down from 
the illumination at the center of the aperture. (It should be noted 
that these descriptive terms, linear tapers, square law tapers, etc. 

“pertain to the field intensity distribution in the aperture and not to 

the power density distribution, which is proportional to the square 

of the field intensity.) 


METHOD OF SUPERPOSITION — 


NORMALIZED APERTURE AMPLITUDE A(r.) — 


0 0.8 0.6 0.4 0.2. 0.0 


NORMALIZED APERTURE RADIUS < 


Figure 3 


Approximation of aperture distribution by uniformly illuminated sub-apertures 
with N = 10. 
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While the on-axis fields have been determined in previous 
literature on this subject for many of the distributions contained 
herein. (see, for example, reference (e)), these fields have been 
recalculated, along with a number of new distributions, using the 
methods of this report. This fulfills several functions. First,. 
the resulting fields of these common distributions are presented, 
for, a number of edge-to-center illumination ratios, in one paper, 
thus providing a convenient comparison of the resulting fields; 
second, it provides a basis for comparing the results of these 
methods with those obtained in previous literature, verifying that 
the results are identical for every distribution for which the 
on-axis field has been previously calculated; and finally, these 
distributions simulate a large number of the possible tapers which 
may occur in antennas, and by choosing the appropriate distribution 
in this report a good approximation of the on-axis field for the 
aperture distribution of interest may be obtained. 


. Linear tapers and square law tapers are two of the most 
commonly used aperture distributions in which the amplitude is 
increased from a low value at the aperture edge to a maximum at the 
center of the aperture. Linear amplitude distributions are given 
by the equationa(r) =1- (1 -a)r where a(r) is the normalized 
amplitude taper, r is the normalized radial. distance from the © 
aperture center and ais the ratio of the amplitude at the aperture 
edge to that at the center. Using similar notation, square law 
distributions may be ‘represented by Q(r) =1- (1 -.a)r*. 


The two families of amplitude tapers, (1 - r=) and (cosine)” 
cover a wide range of possible aperture distributions in which the 
amplitude is gradually decreased from a maximum at the center to a 
somewhat lower value at the aperture edge. The rate of decrease or 
slope may be increased by increasing the value of the parameter n. 
The. on-axis fields have been calculated herein for n = 1/2, 1, 3/2, 
2, 5, and 4 at four different edge-to-center illumination ratios for 
the (cosine)" tapers, and for n =1, 2, 3, and 4 with one illumina- 
tion ratio (o db) in the case of the tapers @(r) = (1 r7)". The 
(cosine)" tapers are expressed by @(r) = (1 - @) (cos S ry>, 

Taylor aperture distributions for circular aperture antennas 
offer an appreciable improvement in beam width and sidelobe levels 
over conventional distributions, and are thus likely to be utilized 
in antenna design. One Taylor distribution is included to provide 
an example of the form of the on-axis field resulting from distribu- 
tions of this type. This Taylor distribution, fram reference Ce). 
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produces a sidelobe level 25 db down from the intensity on the axis 
and has a transition integer of i = 8. Since this is not a mono- 
tonic distribution, the method of finite differences was used to 
calculate the on-axis field. 


Reflector antennas frequently incorporate a primary feed at 
the aperture center which blocks out the center portion of the 
aperture illumination. To determine the effect blocking will 
have upon the field, 5 and wo db cosine tapers with zero illumina- 
tion for r < 0.1 were included. For these distributions also, the 
method of finite differences was utilized. 


Four irregular distributions which could be approximately 


simulated by 5 db cosine?/*, db cosine?/*, 5 db cosine* and 

o db cosine* are included to emphasize the usefulness of the methods 
of this report in determining the fields of irregular distributions 
and to determine the extent of the deviation occurring in the 
on-axis fields when these approximations are made. 


RESULTS 


The on-axis fields were determined by computer for the aperture 
distributions discussed in the previous section and the results are 
presented in the figures of penecae A in a manner believed to Peat. 
facilitate comparison: 


The on-axis fields for the aperture distributions (linear, 
square law and cosine™ tapers) which have four edge-to-center 
illumination ratios are presented with the four resultant fields for 
each taper on the same graph (Figures 5 through 13). In addition, 
Figures 14 through 22 present the (cosine)" tapers with the edge-to- 
center ratio held constant (at each of the values 5 db, 10 db, 20 db 
and o db down from that at the aperture center) and the taper varied 
by letting n take on values 1/2, 1, 3/2, 2, 3 and 4. Since for each 
edge-to~center ratio there exist six distribution functions (one for 
each value of n), the fields for each db level are presented. in two 
figures for clarity with n = 1/2, 3/2 and 3 on one graph and n =1, 
2’and 4 on another. From these arrangements both the effect of 
varying the edge-to-center illumination ratio for a given distribu- 
tion function, and the effect of varying the distribution function 
for a given edge-to-center illumination ratio may be observed. 


LO 
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The on-axis fielas resulting from the (1 - r2)P- ae (with 
n=1, 2, 3 and 4), the Taylor distribution and the irregular 
tapers are presented in Figures 22 through 27. All four fields 
from the (1 - r?)" tapers are presented on the same graph. Each 
of the fields from the irregular tapers is plotted with the regu» 
lar taper that is being approximated. 


\ 


The effect, on the oneaxis field, of blocking due to a 
central feed horn was found to be so slight that it would not be . 
observable if the fields were plotted. Thus these tapers were not 
presented in the figures of Appendix A; however, the computer tabu-« 
lations for these distributions are ineluded in Appendix B, listing 
W p/We at various distances x. 


All of the power density fields in the figures of Appendix A 
are plotted on a logarithmic scale using the dimensionless’ coordi« 
nates W sits as ordinates aut x E 3] as abscissae. The correspond 
ing aperture distributions are included in the upper right hand — 
corner oF eaek figure. 


ERROR ANALYSIS 


To obtain an estimate of the error resulting from the use of 
‘these methods (with N = 50), the exact solutions of the Fresnel 
equation for the 5. db and » db square law and (cosine)* tapers were 
computed. These represent the two extremes of deviation, with the 

5 db square law taper deviating least fram uniform illumination and 
the o db (cosine )* having the greatest deviation from uniform illumi- 
nation. It is assumed that the error ocmurring in the calculations of 
the on-axis fields for other tapers will not exceed that of these 
extreme tapers. Both the 5 db and the o db illumination ratios were 
calculated for each of the two tapers to determine the effect of 
edge-to-center illumination ratios upon the error. 


‘The accuracy of W,/W, in the.tables of Appendix B is limited to 
three significant figures “since the computer input data (the e,s) 
. were only given to three decimal places. Comparison of the exact 
‘values of W p/W, with those given in Appendix B at various values 
‘of x fram 0. 1 to 100, showed that for the majority of cases there 
was no error occurring in the first three significant figures. In 
no instance did the superposition method give an error in excess of 
1.0% and only in the most extreme distribution (o db cosine*) did 
‘the method of finite differences have an error greater than this. 
This amounted to 1 = 3% and can be expected to occur when this 
method is used for the more extreme distributions ( db cosine®, 
cosine* and (1 = r7)*). In both methods the error is greatest at 
low values of x. 


amie 


NWL REPORT NO. 1893 


It is also of interest to observe the effect that varying 
the value of N has upon the accuracy of the calculations. While. 
the method of superposition gives the exact solution of the 
Fresnel equation for a uniformly illuminated aperture, there is 
@ slight amount of deviation when the method of finite differences 
is used. The oneaxis field for a uniformly illuminated aperture 
was calculated using values of N = 10, 25, 50, 75, 100, 150, 200, 
500, and 1000. The resulting values of W/W, were ccupaned with 
the exact solution and the percentage of error was plotted as a 
function of N for various values of x in Figure 4. 


As can be seen from the graph, a linear relationship exists 
(on a logarithmic scale) between the error and N, and if it is 
assumed that other distributions will have a similar error for 
the same value of N, Figure 4 may be used to choose the appropriate 
value of N for whatever accuracy is desired. 


The curves of Figure 4 are applicable only to the method of 
finite differences, since no error occurs when the method of super-= 
position is used for uniform illumination. 


CONCLUSIONS 


Examination of the resulting fields reveals a number of 
significant effects of non-uniform aperture distributions upon the 
on-axis field. These may be summarized as follows: 


1. The deep nulls which occur in the Fresnel field of a 
uniformly illuminated aperture are filled in considerably when the 
distribution is tapered. This is most pronounced for the distribu- 
tions which deviate most from uniform illumination. 


2. There is a decrease in power density in the far field (a 
shifting of the curve to the left) although the slope of this far 
field (inverse square law) curve remains the same. Upon closer 
analysis it is found that the values of W p/Wa in the far field are 


precisely the values predicted by the far field formula 


Mien ON” aiaee. UG = Ue Me 2) a az] 
Wa 40 Re ic A (r) ar 


Gy and G, |- ee 4x] being the gain of the aperture for tapered and 


uniform illumination, respectively. 
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Figure 4 


Error as a function of the number of increments N for uniform illumination 


using the method of finite differences. 
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3. The envelope of the maxima in the Fresnel field of a uni-~ | 
formly illuminated aperture is a constant value, i.e, W/W, = 4.0, 
while for a tapered aperture distribution the peaks tend to either 
increase or decrease as x decreases depending upon the particular 
distribution. From an analysis of the resulting curves and fram © 
theoretical considerations, a formula for this envelope may be 
found, namely that as x approaches zero the maxima approach .the 
value 


. (2) 2 Me (1 a | ee Cr a ee Ls : 
av max Wa , 


2s. a (r) rar 


For tapers which approach W, very rapidly fram the value at 
the aperture edge (square law and (cosine)+/* tapers), the maxima 
decrease in field strength as x decreases, while for most other 
tapers the field strengths of the maxima increase as x decreases. 


4. As the edge-to-center ratio is decreased, the values of x 
at which the near field maxima and minima occur decreases, 1l.e., 
there is a slight shift of the peaks and nulls to the left. 


5. The edge-to-center illumination ratio has more effect upon 
the on-axis field than the particular type of taper used. 


6. "Blocking" of the center portion of an aperture (from 
center out to r $0.1 4) has only a small effect upon the on-axis 
field. The far field is decreased by about 1%, while the near field 

peaks are increased on the order of 54%. 


SUMMARY 


1. Two equations which lend themselves well to computer solu. 
tions have been developed, predicting the on-axis power density for 
any arbitrary aperture distribution. 


2. These equations have been programmed for computer solution 
and are readily available for use as a tool in the study of radar 
fields. , , 
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3. The on-axis flelds were calculated for a number of aper- 
ture distributions which were felt to approximate many of the dis- 
tributions occurring in practice. By choosing the distribution in 
this report which best simulates an actual distribution under BEuey 
a good approximation of the on-axis field may be obtained. 


' 4, Analysis of the resulting on-axis fields reveals a number 
of significant effects which tapered distributions have upon the 
field as compared to uniform distribution. These are: 


a. The nulls occurring in the near field are partially or 
completely filled in. 


b. The far field gain is reduced. 


e. The envelope of the Fresnel field peaks is no longer 
constant but increases or decreases as x decreases. 


d. The values of x at which the Fresnel field peaks 
occur decrease. 


e. The edge-to-center illumination ratio has more effect 
upon the on-axis field than the particular type of taper used. 


£. Blocking of a small center portion of the aperture 
(r = 0.1 a) has little effect upon the on-axis field. 
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